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We present 2D simulations of pair-instability supernovae considering rapid rotation dur¬ 
ing their explosion phases. Recent studies of the Pop III star formation suggested that 
these stars could be born with a mass scale about 100 Mq and with a strong rotation. 
Based on stellar evolution models, these massive Pop III stars might have died as highly 
energetic pair-instability supernovae. We perform 2D calculations to investigate the im¬ 
pact of rotation on pair-instability supernovae. Our results suggest that rotation leads 
to an aspherical explosion due to an anisotropic collapse. If the first stars have a 50% of 
keplerian rotational rate of the oxygen core before their pair-instability explosions, the 
overall 56 Ni production can be significantly reduced by about two orders of magnitude. 
An extreme case of 100% keplerian rotational rate shows an interesting feature of fluid 
instabilities along the equatorial plane caused by non-synchronized and non-isotropic 
ignitions of explosions, so that the shocks run into the in-falling gas and generate the 
Richtmyer-Meshkov instability. 
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1. Introduction 

The first generation of stars so-called, Population III (Pop III) stars were predicted 
to form about several hundred million years after the Big Bang [U [2]. The Pop III 
stars were born in the gravitational-potential wells constructed by the dark matter 
halos of masses about 10 6 solar masses (Mq), that allowed the primordial gas (mass 
content: 76 % hydrogen, 24% helium) to form stars. Because no heavy elements 
(metals) were present at that time, molecular hydrogen served as the most efficient 
coolant, but it could not cool the gas temperature efficiently. Thus, the Jeans mass 
of the Pop III star-forming cloud was expected to be much larger than the one of 
the current star formation regions; consequently the Pop III stars could be more 
massive than the present-day stars. Modern cosmology simulations suggested that 
the mass scale of the Pop III stars would fall between tens to hundreds of Mq. 
The exact mass scale is still under debate. Without initial metals such as carbon 
and oxygen, the Pop III stars needed to reach higher core temperatures to burn 
hydrogen effectively |3J, that also resulted in higher surface temperatures. So these 
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stars could produce large amounts of UV photons and extensively ionized primordial 
hydrogen, and helium inside the inter-galactic medium (IGM) leading to the cosmic 
reionization After their short life time of several million years, many 

massive Pop III stars might have died as supernovae and dispersed the first metal 
to the primordial IGM. Such chemical enrichment might have changed the later star 
formation mode from Pop III to Pop II and profoundly reshaped the simple early 
universe into a state of ever-increasing complexity. Thus, understanding the first 
stars and their supernovae has become one of the frontiers of modern cosmology. 

Fate of a massive star strongly correlate to its final mass before its death. Be¬ 
cause there is no initial metal inside the envelope of Pop III stars, the mass loss 
rates due to the metal line-driven wind can be strongly suppressed. Hence, it is rea¬ 
sonable to assume that Pop III stars losses very little or no mass during their stellar 
evolution. Based on stellar physics mm, the Pop III stars with initial masses of 
10 — 80 Mq eventually forge an iron core about one solar mass. When the iron core 
mass exceeds its Chandrasekhar mass m, the degenerate pressure of electrons is 
no longer able to support its own gravity. The iron core catastrophically collapses 
into a black hole or neutron star. The energy released from the gravity allows an 
energetic explosion so-called a core-collapse supernova. However, such explosions 
may fail and the entire star may collapses directly into a black hole. The nature of 
core-collapse supernovae is complicated by several hurdles such as neutrino physics, 
multi-scale, and multi-dimension nu [in nn nn m]. The explosion mechanisms re¬ 
main subject to investigation. If Pop III stars are more massive than 80Mq. after 
the central carbon burning, their cores encounter the electron/positron production 
instability, the core reaches sufficiently high temperatures (~ 10 9 K) and low den¬ 
sities (~ 10 6 gem -3 ) to favor the creation of electron-positron pairs. It reduces the 
adiabatic index 7 below 4/3 and causes a dynamical instability of the core. The core 
is now at quasi hydrodynamical equilibrium and its temperatures begin to oscillate 
for a period of its dynamic time scale of several hundred seconds. The variations 
in temperature do not cause a significant impact to its evolution and the star later 
dies as a core-collapse supernova. However, if the stars are more than 100 Mq, the 
oscillation of temperatures becomes very violent. Several large temperature spikes 
produce shocks with energy of 10 49 — 10 5 °erg. These shocks can not disperse the 
entire star, but eject its outer envelopes instead. Due to different amount of mass 
and energetics, these ejected shells are likely to run into each other. The colliding 
shells effectively convert the kinetic energy of ejecta into the internal energy of gas. 
The collision sites usually happen at optically thin regions and thermal radiation 
escape immediately and produce very luminous optical transits known as pulsa- 
tional pair-instability supernovae mm- Once the mass of stars reaches between 
150 and 260 Mq. At this point, pair-instabilities become so violent and trigger a 
runaway collapse of core. The core temperature and density rise swiftly and ignite 
the explosive oxygen/silicon burning. The energy released from the burning converts 
the contraction of the core into an energetic explosion that completely disrupts the 
entire star. This thermonuclear explosion, known as a pair-instability supernova 
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(PSN), that has explosion energies up to 10 53 erg and yield up to 50 M 0 of 56 Ni 
[m [201 EH |22l EH ISl 125] . Most of Stars are 300 M 0 are believed to simply die as a 
black hole. A new type of supernova caused by the general-relativity instability in 
supermassive stars at ~ 55,000 M 0 , which may have formed in the early universe, 
has now also been found [26]. 

Because the Pop III stars were predicted to very massive, many of them tended 
to die as PSNe. The ID PSN models have been extensively studied by several 
authors [m [in uni The extensive yields, explosion energies, and radiation 
properties of PSNe are calculated in their models. Unlike core-collapse supernovae, 
the thermonuclear explosions of PSNe are not very sensitive to the dimensionality 
of simulations. However, simulating the mixing of the fluid instabilities in PSNe still 
requires multidimensional simulations that recently become accessible and studied 
by [28] HU [221 281 • The simulations of [30] found the mixing of PSN from the both 
burning and hydrodynamics instabilities. In the red super-giant progenitors, the 
mixing driven by the Rayleigh-Taylor instabilities can be very significant. 

One of the important physical quantities in the stellar evolution is rotation. In 
the real world, the process of star formation usually endows a certain amount of 
angular momentum to the star and makes it rotate. Due to the technical difficulty 
of the study, the rotational rates of massive stars are still poorly understood in 
both theoretical and observational studies. Rotation can affect the stellar evolution 
as well as the resulting supernova explosions. The nucleosynthesis inside the stars 
changes due to the rotational mixing, so that the chemical composition of the star 
changes. If metal from the inner part is mixed out to the envelope, the stellar wind 
can be enhanced due to metal-driven lines. Strong winds reduce the stellar mass 
dramatically and affect the fate of the stars. Studies of the stellar evolution of very 
massive stars with rotation have been performed by [20] 31, 32] [33] 134]. Recent 
results from [35] [36j show that rotation can lower the mass criterion for PSNe 
progenitors because mixing facilitates helium burning, resulting in a more massive 
oxygen core. Recent simulations m suggested that the Pop III stars could rotate 
very fast, up to 50% of keplerian rate at the surface. If such a high rotation does 
exist, it can affect the evolution of the stars and their supernovae. So we perform 2D 
simulations of PSNe with rotation and investigate how the strong rotation impacts 
the explosive burning. The purpose of this paper is to explore the energetics and 
nucleosynthesis of PSNe during the explosion phase. 

The structure of this paper is as follows: we first describe the numerical ap¬ 
proaches and setup of our simulations in § [2] Then we present the results in § [3] and 
conclude our findings in § [4j 

2. Methodology and Problem Setup 

The 3D stellar evolution models followed from the main sequence star to supernovae 
are still unavailable due to the limitations of current computational power. Alterna¬ 
tively, we start simulating the stars with KEPLER [SHIES], a ID spherically-symmetric 
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Lagrangian code, and follow the evolution of stars up to 20 — 100 seconds before 
the supernova explosions. Then we map the resulting ID stellar evolution models 
from KEPLER onto 2D grids of CASTRO as the initial conditions. We then apply our 
differential rotation scheme to CASTRO and follow the simulations until the explosive 
burning ceases when the shock has successfully launched. This setup is designed to 
model the most critical phases of the supernova explosion in 2D with practical com¬ 
putational resources. In this section, we introduce our progenitor models, problem 
setup, and numerical methods. 


2.1. Progenitor Models 

We use progenitor Pop III stars of 150 M 0 , 200 M 0 , and 250 Mq with very little 
overshooting at the late-time evolution. These stars eventually become blue super¬ 
giants which have smaller radii than those of red supergiants. The ID stellar models 
are evolved until most of the explosive burning is about to happen. Physical prop¬ 
erties of these stars are listed in Table 1. The resulting ID profiles are then mapped 
onto 2D cylindrical grids of CASTRO using a new mapping algorithm [40j. This ini¬ 
tialization scheme conservatively maps the physical quantities such as mass and 
energy during mapping from ID to multi-D. The initial perturbations are seeded in 
the form of velocities based on the stellar convection physics. 


Table 1. ID Progenitor models 


Name 

M* 

[Mq] 

AT He 

[M @ ] 

pc 

[10 6 g cm -3 ] 

T c 

[10 9 K] 

R 

[10 13 cm] 

B150 

150 

67 

1.40 

3.25 

16.54 

B200 

200 

95 

1.23 

3.31 

2.86 

B250 

250 

109 

1.11 

3.34 

23.06 


M *: initial stellar mass; Mjj e : helium core mass at collapse; p c : central density at 
collapse; T c : central temperature at collapse; and R: stellar radius. 


2.2. CASTRO 

We run 2D simulations using CASTRO j4T[[42], a massively parallel, multidimensional 
Eulerian, adaptive mesh refinement (AMR), hydrodynamics code for astrophysical 
applications. CASTRO was originally developed at the Lawrence Berkeley Lab, and it 
is designed to run effectively on supercomputers over 10,000 CPUs. CASTRO provides 
a powerful platform for simulating hydrodynamics and gravity for astrophysical gas 
dynamics. Modeling thermonuclear supernovae requires calculating the energy gen¬ 
eration rates from nuclear burning, which occurs over an extensive range of temper¬ 
atures, densities, and compositions of isotopes. We have implemented the APPROX 
19-isotope reaction network [35, ]33j into CASTRO. It is capable of efficiently calcu¬ 
lating accurate energy generation rates for major nuclear reactions from hydrogen 
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to silicon burning. In CASTRO, we use a realistic equation of state (EOS) [Uj for 
stellar matter. This EOS considers the (non)degenerate and (non)relativistic elec¬ 
trons, electron-positron pair production, as well as ideal gas with radiation. It is a 
tabular EOS that reads in p, T, and X-, of gas and yields its derivative thermody¬ 
namics quantities. CASTRO offers different types of calculation for gravity, including 
Constant, Poisson, and Monopole. For supernova simulations, spherical symmetry 
is still a good approximation for the mass distribution of gas at the early phase of 
evolution. Using this approximation saves lot of computational time in the gravity 
solver. For 2D or 3D CASTRO simulations, we first calculate a ID radial average 
profile of density, then compute the ID profile of g field and use it to calculate the 
gravity of the multidimensional grid cells. Care is taken to resolve the important 
scales of the explosions such as catching the shock front and mixing driven by fluid 
instabilities. The grid structure for 2D simulations uses base grids, 256 x 512 with 
additional three levels of AMR to resolve a domain of (4x 10 11 ) x (8x 10 11 ) cm 2 . This 
setup yields the size of the finest patch about (7.2 x 10 7 ) 2 cm 2 which fully resolves 
the nuclear burning and the structures of fluid instabilities of rotating PSNe. The 
AMR criteria are set for gradients of density, velocity, and pressure. Finer zones 
are automatically created on the top of over-gradient regions. 2D CASTRO uses a 
cylindrical coordinate r — z. The axis z serves as the rotational axis. Since we sim¬ 
ulate only half of the star in 2D. The lower boundary of r uses reflect conditions to 
prohibit fluid from entering; the other three boundaries use outflow conditions that 
allow the fluid to freely cross over. We also constructed nested zones to ensure that 
the inner region of the star constantly receives a higher spatial resolution. 


2.3. 2D Rotational Model 

We develop a differential rotational model motivated by the stellar structure. It 
yields non-uniform rotational rates inside different regions of stars. Because the 
oxygen core of the star is a dense and compact object, we assume it behaves like a 
rigid body with a constant rotational rate, w. Outside the oxygen core, the structure 
of the helium and hydrogen envelope is relatively puffy and unlikely to keep the 
constant ui. Instead, we assume a constant specific angular moment extending from 
the edge of oxygen core to its surface and w monotonically decreases from the 
rotational axis, as shown in Figurc[l] To determine the radius of oxygen core and its 
rotational rate, we plot the oxygen abundance as a function of radius in Figure |2 (a) | 
The most abundant oxygen is located at the radius between 10 9 — 10 10 cm, and its 
outer boundaries of the oxygen core is at 10 10 cm. Based on the size and mass of 
the oxygen core, we determine its keplerian rotational rates, u>*. which is about 
0.5 sec” 1 . To examine whether the model is viable, we calculate tOk as a function 
of radius in Figure [2(b)] It shows a nice trend of constant uj in the the oxygen core 
and slowly decrease from the core boundary. Our assumption of oxygen rotating as 
a rigid body seems to be very suitable for the model. 

We therefore assume a constant rotational rate inside the oxygen core and this 
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Fig. 1. Illustration of the differential rotational model, uj is the rotational rate and r is the distance 
to the rotational axis, uj inside the oxygen core (r < r c ) is assumed as a constant. Beyond the 
core (r > r c ), we assume a specific angular j constant and result in uj oc r ~ 2 . 



(a) Oxygen abundance pattern (b) uj & as a function radius 

Fig. 2. (a) The size of the oxygen cores of presupernovae can be identified by their oxygen 

abundance. The sudden drop in the oxygen core marks the size of oxygen core, (b) uj{r) behaves 
like a constant then decreases. The saddle point of the drop is roughly at the edge of oxygen core. 


rate is based on the keplerian rate of the oxygen core. In CASTRO, the rotation is 
calculated by the angular momentum, j which provides a force term in the Euler 
equation for fluids and itself is evolved by advect equations: 


d(pu) 

dt 

d{pj) 

dt 

F c 


-V • (puu) - Vp + pg + F c , 

-V • (puj). 
j 2 

p'-^\ j = wr 2 when r < r c , otherwise j = ur 2 


Cl) 

( 2 ) 


(3) 
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Here p,p, u, g are the density, pressure, velocity vector, and gravitational vector, 
respectively. F c is the centrifugal force calculated by p, the angular momentum 
per unit mass j, and the distance to the rotational axis r z . j is only initialized at 
the beginning of the simulations, then it evolves with fluid elements following the 
equations as above. 


3. Results 

We first present the results from the representative model of the 200 Mq rotat¬ 
ing Pop III star with w = 0.5wfc. After the onset of the 2D CASTRO simulation, 
the core of the star starts to collapse due to its dynamical instabilities caused by 
pair-production instabilities. Similar to the case of non-rotating stars, explosive 
oxygen and silicon burning proceed as the core contracts. After about 25 seconds, 
the energy released from the nuclear burning launches a shock and the core of star 
has successfully exploded. For a non-rotating case, collapse of the core is perfectly 
anisotropic and homogeneous. Thus, the resulting explosion appears to have spher¬ 
ical symmetry. However, for a rotating star, rotation provides a centrifugal force 
to resist the in-falling gas during the collapse. Because the centrifugal force follows 
as F c oc r~w 2 , for the given distance from the center of star, the core receives the 
strongest F c along the equatorial direction and no force in the polar direction. It 
automatically leads to anisotropic compression of the core and reflects in the explo¬ 
sion. The maximum compression happens along the polar direction and yields more 
explosion energy and releases a more powerful shock wave than other directions. 
The core is blew up in an elliptical shape with a long axis along the rotational axis. 
Figure [3] shows the post-explosion of 200 Mq progenitor star, ui = 0 in the left-half 
panel and ui = 0.5wfc in the right-half panel. Both snapshots are taken at the same 
time after the explosions occur. In comparison of both results, the non-rotating 
model shows a stronger shock wave, and demonstrates some mixing at the inner 
part of the oxygen-burning shells. The rotating model, however, shows a relatively 
less energetic explosion due to a weaker compression; the ejected oxygen-burning 
shell is shaped into an ellipse and shows very little mixing. PSNe of non-rotating 
stars usually produce a large amount of radioactive isotope, 56 Ni which is made 
inside PSNe mainly by the explosive 28 Si burning right before the core bounces. 
Because the decay energy from 56 Ni powers the PSN light curves, it becomes one 
of the key isotopes for observations. The amount of 56 Ni production can determine 
whether the PSNe become bright or faint SNe. We show the corresponding 56 Ni 
abundance and gas density in Figure |4j The non-rotating star demonstrate a visual 
56 Ni abundance, which is missing in the rotating model. It suggests that 56 Ni pro¬ 
duction has been significantly suppressed due to rotation. Only about 10 -2 Mq 56 Ni 
is made in the rotational model. However, the non-rotating model has synthesized 
about 6.57 Mq of 56 Ni. This suggests that a strong rotation can significantly affect 
the energetic and 56 Ni production in PSNe . 

What happens if ui = ujk ? This is an extreme rotational rate which can al- 
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Fig. 3. Comparison of rotating and non-rotating models from a 200 Mg Pop III star. The left-half 
sphere is the non-rotating model, and the right-half is the rotating model of u ~ 0.5a.’/,. The gray 
colors show the radial velocities v r , and the colorful contours are the 16 O mass fraction. Both 
snapshots are taken at the same time, about 200 sec after the onset of the explosion. 


most break out of the star’s core. To answer this question, we simulate a 200 Mq 
of to = uik- The anisotropic compression becomes even more stronger than that of 
the u> = 0.5wfc case. The ellipse shape of the ejected core has an larger eccentricity. 
One interesting fluid instability has been found in the carbon-burning shell, which 
is shown in Figure [5j The carbon shell breaks along the equatorial plane and over¬ 
shoots some carbon into the oxygen core. This overshooting develops strong fluid 
instabilities. Because the explosion is anisotropic, the shocks of different directions 
are initialized at different times: the pole comes first, and the equator comes last. 
In this fast-rotating model, the shock from the pole has been sent out, but the gas 
along the equator is still collapsing. Once the shock runs into the collapsing gas, 
the Richtmyer—Meshkov (RM) instability [33] starts to develop and the resulting 
mixing makes carbon penetrate into the oxygen-rich envelope. 
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Fig. 4. Comparison of rotating and non-rotating models from a 200 Mg star. The left-side sphere 
is the non-rotating model, and the right side shows a rotating model of u) ~ 0.50;*,. The gray color 
shows the density, and the colorful contours are the 5,1 Ni mass fraction. Both snapshots are taken 
at the same time, about 200 sec after the onset of the explosion. 


4. Conclusions 

Stellar rotation plays a key role in both stellar evolution and the fate of the very 
massive Pop III stars. The nature of the rotation can be traced back to the phase 
of star formation, when they are endowed with angular momentum. The definite 
rotational rates for the massive Pop III stars are still unknown. Results of these star 
formation suggest that a high rotational rate is 50% of its keplerian rotational rate. 
We present the results of the impact of rotation on a 200 M 0 star. The 50% keplerian 
rotational rate case demonstrates the onset of an anisotropic explosion and strongly 
suppresses 56 Ni. If this does apply for PSN progenitors, 56 Ni production shown in 
previous PSN models may change and the corresponding PSNe luminosity will be 
attenuated. An extreme case of a 100% keplerian rate shows an interesting feature of 
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Fig. 5. Hot-color shows the mass fraction of 12 C. Strong fluid instabilities occur during the 
explosion along the equatorial plane shown in the red box of the left panel; its close-up appears in 
the right panel. 


overshooting along the equatorial plane caused by non-synchronized explosions that 
send shocks into the infilling gas. Finally, we conclude that PSNe of fast rotating 
stars would look much fainter than we originally expect because of the limited 56 Ni 
production. 
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